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Abstract. We examine the multiple Hamiltonian structure and construct a 
symplectic realization of the Volterra model. We rediscover the hierarchy of 
invariants, Poisson brackets and master symmetries via the use of a recursion 
operator. The rational Volterra bracket is obtained using a negative recursion 
operator. 



1. Introduction 
The Volterra model, also known as KM system is defined by 

(1) iii = Ui(u i+ i - Ui-i) i = l,2,...,n, 

where uq = =0. It was studied originally by Volterra in [T7] to describe 

population evolution in a hierarchical system of competing species. It was first 
solved by Kac and van-Moerbeke in [13], using a discrete version of inverse scat- 
tering due to Flaschka [TUj- In [TS] Moser gave a solution of the system using the 
method of continued fractions and in the process he constructed action-angle co- 
ordinates. Equations ([T]) can be considered as a finite-dimensional approximation 
of the Korteweg-de Vries (KdV) equation. They also appears in the discretization 
of conformal field theory; the Poisson bracket for this system can be thought as a 
lattice generalization of the Virasoro algebra [5] . The variables itj are an interme- 
diate step in the construction of the action-angle variables for the Liouville model 
on the lattice. This system has a close connection with the Toda lattice, 

a,: = <Zi(&i+i - hi) % = 1, . . . ,n - 1 

k = 2(a? - af_j i = l,...,n. 

In fact, a transformation of Hcnon connects the two systems: 

a-i = -^u 2 iU 2 i-i i=l,...,n-l 

h = ^("2i-i + u 2l -2) i = l,...,n. 

We note that the number of variables for the Toda lattice is odd and therefore 
we restrict our attention to the Volterra system with an odd number of variables. 
The Volterra system is usually associated with a simple Lie algebra of type A n . 
Bogoyavlensky generalized this system for each simple Lie algebra and showed that 
the corresponding systems are also integrable. See [TJ [2] for more details. The 
relation between Volterra and Toda systems is also examined in [5]. 

The Hamiltonian description of system ([T]) can be found in [7] and [3]. We will 
follow [3] and use the Lax pair of that reference. The Lax pair is given by 

L=[B,L], 

where 
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This is an example of an isospectral deformation; the entries of L vary over time 
but the eigenvalues remain constant. It follows that the functions Hi = iTrl? are 
constants of motion. We note that 



Hi 



corresponds to the total momentum and 

n 



Hn 



i=l 



n-1 

! 22 UiU i+1 

i=l 



is the Hamiltonian. 

Following [3] we define the following quadratic Poisson bracket, 

{Ui, U i+ i} = UiU i+ l, 

and all other brackets equal to zero. We denote this bracket by ir 2 . For this bracket 
detL is a Casimir and the eigenvalues of L are in involution. Of course, the functions 
Hi are also in involution. Taking the function y]" Ui as the Hamiltonian we obtain 
equations HJ . This bracket can be realized from the second Poisson bracket of the 
Toda lattice by setting the momentum variables equal to zero [7] . 

In [5] one also finds a cubic Poisson bracket which corresponds to the second 
KdV bracket in the continuum limit. It is defined by the formulas, 

{u tl U. l+1 } = UiU i+ l(u t + Uj+l) 
{u l: U l+2 } = UiU i+ lU l+2 , 

all other brackets are zero. We denote this bracket by ir 3 . In this bracket we still 
have involution of invariants. We also have Lenard type relations of the form 

7T 3 VH L =1T 2 VH 1+1 . 

In [3] appears a bracket that is homogeneous of degree one, a rational bracket 
constructed using a master symmetry. This bracket, denoted by m, has TrL as 
Casimir and the Hamiltonian is H 2 = |TrL 2 . The definition of the bracket is the 
following. We define the master symmetry F_i to be 
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where the /j are determined recursively as follows, 



fl — — 1) /2i — /2i-li / 2 i-l — —fli-1 — 1- 

«2i-l 

Taking the Lie derivative of tt 2 in the direction of Y-\ we obtain 7Ti, a Poisson 
bracket that is homogeneous of degree 1. For n = 5, 7Ti takes the form: 



W,it 2 } 


= u 2 


{Ui,U 3 } = 


-u 2 




U2U4 

"3 

U2U4 
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{mi,w 5 } 
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= u 2 


{u 2 ,u 4 } = 


U2U4 

U3 


{u 2 ,u 5 } = 




{u 3 ,u 4 } 


= U4 


{u 3 ,u 5 } = 


—U4 


{u 4 ,u 5 } = 


U4. 





In this paper we rediscover this bracket using a recursion operator. The higher 
Poisson brackets are constructed using a sequence of master symmetries Yi, i — 
— 1,0,1, ... . We define Yq to be the Euler vector field 



n d 
i— 1 

The explicit formula for Y\ is 

2 — 1 

where 

Ui = (i + l)uiu i+ i + uf + (2 - i)ui-iUi . 

It is easily checked that the bracket ir 2 is obtained from 7Ti by taking the Lie 
derivative in the direction of Y\ . Similarly, the Lie derivative of ~k 2 in the direction 
of Y\ gives 7r 3 . 

The brackets m, -k 2 and ir^ are just the beginning of an infinite hierarchy con- 
structed in [3] using master symmetries. We quote the result: 

Theorem 1. There exists a sequence of Poisson tensors itj and a sequence of 
master symmetries Yj such that: 

i) iTj are all Poisson. 

ii) The functions Hi are in involution with respect to all of the ttj. 
in) Y % {H ) = {i+j)H i+3 . 

iv) L Yi TTj = (j -i - 2)TT l+j . 

v) [Y h Y J ] = (j-i)Y i+j . 

vi) 7TjV-ffj = 7Tj_iVi?i+i, where itj denotes the Poisson matrix of the tensor Ttj. 

In this paper we prove the results of Theorem 1 using a different approach. 
Namely, we construct a recursion operator in a symplectic space, define all master 
symmetries, invariants and Poisson brackets using results of Magri and Oevel and 
then project to the space of u variables. 
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2. Master Symmetries and a Theorem of Oevel 

We recall the definition and basic properties of master symmetries following 
Fuchssteiner [12]. Consider a differential equation on a manifold M defined by a 
vector field X - We are mostly interested in the case where x is a Hamiltonian vector 
field. A vector field Z is a symmetry of the equation if 

[z, x ] = o. 

A vector field Z is called a master symmetry if 

[[z,x],x] = o, 

but 

[z,x]?o- 

Master symmetries were first introduced by Fokas and Fuchssteiner in [11] in con- 
nection with the Benjamin-Ono Equation. 

A bi-Hamiltonian system is defined by specifying two Hamiltonian functions Hi , 
i?2 and two Poisson tensors tt± and 7T2, that give rise to the same Hamiltonian 
equations. Namely, 7riVi?2 = K2VH1. The notion of bi-Hamiltonian structures is 
due to Magri |14j . Suppose that we have a bi-Hamiltonian system defined by the 
Poisson tensors 7Ti, 1^2 and the Hamiltonians Hi, Hi- Assume that 7ri is symplectic. 
We define the recursion operator 1Z — 7r 2 7rf , the higher flows 

Xl = TV- 1 xi , 
and the higher order Poisson tensors 

m = tv^tti . 

For a non-degenerate bi-Hamiltonian system, master symmetries can be generated 
using a method due to Oevel [IB] . 

Theorem 2. Suppose that Xq is a conformal symmetry for both m, ^2 and Hi, 
i.e. for some scalars X, \x, and v we have 

Cx ^i = Xir 1 , Cx K2 = ^2, Cx Hi = vH\. 
Then the vector fields Xi = TZ 1 Xq are master symmetries and we have, 
(a) Cx.Hj = (v+( 3 -l + i){fx - X))H l+J 
(6) C x ,ttj = (m+ (J ~ i ~ 2)(/x - X))n l+j 
(c) [X i ,X j ] = (»-X)(j-i)X i+j . 

3. Symplectic realization 
We define the following transformation from R 2 ™ to R 2n_1 , 
U21-1 = -e Pl i = 1, . . . ,n , 

(3) U2 l = e q ^ 1 - q * i = l,...,n-l. 
The Hamiltonian in (q,p) coordinates is given by 

n n— 1 

(4) hi = -J2 ePz + J2 e "* +1 ~ qz - 

i=l i=l 

It is straightforward to check that Hamilton's equations for (£[]) correspond in the 
u— space to the KM-system ([T]) via the mapping Q. The symplectic bracket in 
(q, p) coordinates corresponds to the quadratic bracket ir 2 ■ For this reason we will 
denote the standard symplectic bracket in R 2 ™ by J2. Our purpose is to define a 
bracket J3 in R 2 ™ which is mapped to 7T3 under the transformation ([3]). The idea 
of the construction is to lift the master symmetry Yi from the it—space up to the 
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{liP)~ space and obtain a vector field which we denote by X\. The new bracket 
J3 will be defined as the Lie derivative of J 2 in the direction of Xi . One possible 
definition for X\ is the following: 

n £) n £) 

*-£*£ + £< 

where, 

i-l i-1 
Aj = -e Pl - ^ & + (1 - 2 z)e p * + ^ e «*+ 1 -«* i = 1, 2, . . . , n , 

B l =2i e qi + 1 - qi - e p ' + (3 - 2 ije 9 * -9 *- 1 i = 1, 2, . . . , n . 

We note that in the summations if an index is not defined then we ignore that 
whole term. 

Taking the Lie derivative of the symplectic bracket J 2 in the direction of X\ we 
obtain the Poisson bracket J3, 

(5) { qi ,qj} = e p * l<j<t-l<n-l 

= + e 9 * - ^- 1 i = l, ...,n 

{q l ,Pj\ = e^-v- 1 - e ^ +1 ~ q i 1 < j < i - 1 

{ Pi ,p i+1 } = e^~ qi i = l,...,n-l 

The Jacobi identity is straightforward to check. There are four cases (three p, three 
q, two p one q, and two q one p). Two of the cases are trivial and the other two 
can be broken-up to at most five subcases. 

Let J 2 be the symplectic bracket with Poisson matrix 





-/ o y 

where I is the n x n identity matrix. The bracket J 2 is mapped precisely to 
the bracket 7T2 under transformation ([3]), and J3 corresponds to 773. We define 
a recursion operator as follows: 

K = JzJz 1 - 

This operator raises degrees and we therefore call it the positive Volterra operator. 
In (q,p) coordinates, the symbol Xi is a shorthand for Xh t - It is generated, as usual, 
by 



Xi = K i - 1 X i. 



For example, 



1 n n— 1 n — 1 

i—l i— 1 z— 1 

Note that h 2 corresponds under mapping (|3|) to a constant multiple of H 2 
^Tr (L) 2 . In a similar fashion we obtain the higher order Poisson tensors 

J,=ft ,; - 2 J 2 £ = 3,4,.... 
We finally define the conformal symmetry 



n n a 
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The Poisson tensors J2, J3 and the functions define a bi-Hamiltonian pair. 

Namely, J2V/12 — J3V/11. We note that J3 is automatically compatible with J2 
since it is constructed using a master symmetry (see [4] p. 5518). It is straightfor- 
ward to verify that 

C'XqJi = 0, £-x Q <h = J3, Cx {h\) = h\. 

Consequently, Xq is a conformal symmetry for J2, J3 and hi. The constants ap- 
pearing in Oevel's Theorem are A = 0, ji = 1 and u = 1. Therefore, we end up with 
the following deformation relations: 

[Xi, hj] = (i+j)h i+j 
L x,Jj = (j — i — 2)J i+ j 

[X i ,X j ] = (j-i)X l v . 

Projecting to the u— space under mapping (J3]) we obtain relations (iii)-(v) of The- 
orem Q] Statements (i) and (ii) of Theorem Q] follow easily from properties of the 
recursion operator. 



4. The negative Volterra hierarchy 

In this section we describe how the first bracket Tti is obtained via the use of the 
negative operator. The negative operator was introduced in [5] in connection with 
the Toda lattice. We define J\ as follows: 

Ji = WJ 2 , where M = J2 J3" 1 . 

We then project the J\ bracket to the u— space using transformation ([3]) to obtain 
the bracket m. We illustrate in detail the case n = 5. 

We consider the Volterra model in R 6 with coordinates (qi, <73,Pi,P2,P3)- 
Transformation ^ is given by, 



(6) ui = -e Pl , u 3 = ~e P2 , u 5 = -e P3 , u 2 = e q2 ~ qi , u 4 = e 93 " 92 . 

j2= (-l 0)' 

where I3 is the 3x3 identity matrix, and J3 is the Poisson matrix §E§ , 



/ -e Pl -e Pl -e Pl 

e Pl -e P2 - e Q2~qi _ e P2 + e 92-9l 

e Pl e P2 _ e 92-9l e 92-9l _ e9 3-92 _gP3 +e 93 

e Pi e 92-gi e 92-9i e q2 ' qi 

e P2 — e q2 ~ qi —e q2 ~ qi + e q3 ~ q2 — e Q2-qi q e <?3-<72 



J3 = 



\ 

One can find the matrix Ji, 



=93-92 
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(Ji)i,2 = -^e Pl (e P3 - e«-« 2 ) (J^ = ^(e P2 - e ft-») 



(Ji)i,4 = -e^[( e p3 - e*" 92 ) - e P2 e P3 ] (Jx)^ = £e»*e«>-® 

(Ji)i,e = --e^e®-® (J!) 2 ,3 = 

(Ji) 2 ,4 = ie^(e P3 - e «-« 2 ) (J!) 2 , s = ie^(e« _ e P 3 ) 

(Ji) 3)4 = -i e ®-« l e*»-» (Ji) 3 ,6 = £e»*e«»-« 

(^1)3,6 = ~e^e™ (J x )4,5 = -£e»e®-« 

(Ji) 4 , 6 = ie 92 -"^ 93 -® (JOs.e = -ie^e"-* 

where £> = e Pl e P2 e P3 . We note that Z? corresponds in the u— space to the square 
root of det(L). The projection of J\ to the u— space under transformation © is 
precisely the bracket m given in equation ((2]), e.g. 

{ux.ua} = {-eP\e<n-«i} = -e^e^-^({ Pl ,q 2 } - {puQi}) 

= - fp^Xepl [e q2 - qi (-e P3 + e^'^) + e q - 2 - qi e P3 - e q2 - qi e q3 ~ q2 - e P3 e P2 } 

e q 

Using the recursion operator J\f we can construct the negative Volterra hierarchy, 
i.e. J,_i = AfJi, i = l,0, —1, —2, — Using the same method of proof as in [5] one 
can easily show that the conclusions of Theorem [T] hold for any integer value of the 
index. For example, for i = 1 we obtain a Poisson bracket Jo, which projected to 
the u— space gives a rational Poisson bracket of degree zero, tto. In the case of the 
Volterra model in R 4 one can find that ttq is given by, 

r , "2(^2 +W3) r i u 2 (u 1 +u 2 + u 3 ) 

{ui,u 2 } = {u 3 ,wi}- 



{U2,U 3 } 



U1U3 U1M3 
1t 2 (?i 2 + Ml) 



5. Conclusions 



This paper contains three main ingredients. The first, consists of the odd- 
dimensional space of the Volterra model together with its multiple Hamiltonian 
structures. These results are not new but they are derived in this paper using an 
entirely new approach. The quadratic and cubic brackets 7r 2 and 7T3 are contained 
implicitly in the book or Fadeev and Takhtajan [7j. The rational linear bracket w± 
and the rest of the hierarchy were first computed in [3] using master symmetries. 

The second part is a realization of the model in a symplectic space. We define 
a Hamiltonian system in (q,p) coordinates, we compute master symmetries and 
a second Poisson bracket which defines a bi-Hamiltonian pair. We then use the 
recursion operator to produce the infinite hierarchy. All the results here are new. 

The third part is the mapping which connects the two spaces and the two sys- 
tems. It is a mapping from an even 2n dimensional space to an odd 2n — 1 di- 
mensional space. This symplectic realization is also new. We have to mention that 
there is another symplectic realization of the model which goes back to Volterra. 
However, the map is from a An — 2 to a 2n — 1 space, see e.g. j5J. Due to the big 
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difference in dimension the resuits of the present paper will be difficult to duplicate 
for that particular realization. 
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